library(foreign)
library(readstata13)

d <- read.csv("Figure_3.csv", header = TRUE)
grass <- read.dta13("Barber_Bolton_Thrower_Replication_Data_LSQ.dta")

plot(d$polar, d$divided, 
	main = "Marginal Effect of Divided Government at Different\nLevels of Polarization",
	xlab = "Polarization", ylab = "Effect of Divided Government",
	xlim = c(0.46, 3.1),
	ylim = c(min(d$lo90), max(d$hi90)), type = "l", lwd = 3)
polygon(c(d$polar, rev(d$polar)), 
	c(d$lo90, rev(d$hi90)), col = rgb(0.5, 0.5, 0.5, 0.5), border = F)
lines(d$polar, d$divided, lwd = 3)
abline(h = 0)
rug(grass$house_polarization)
